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Typically options with a path dependent payoff, such as Target Accumulation 
Redemption Note (TARN), are evaluated by a Monte Carlo method. This paper 
describes a finite difference scheme for pricing a TARN option. Key steps in the 
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\C ' proposed scheme involve tracking of multiple one-dimensional finite difference 

r-^ I solutions, application of jump conditions at each cash flow exchange date, and 

tH- I a cubic spline interpolation of results after each jump. Since a finite difference 

scheme for TARN has significantly different features from a typical finite dif- 



ference scheme for options with a path independent payoff, we give a step by 
step description on the implementation of the scheme, which is not available in 
^ I the literature. The advantages of the proposed finite difference scheme over the 

j^ , Monte Carlo method are illustrated by examples with three different knockout 

types. In the case of constant or time dependent volatility models (where Monte 
Carlo requires simulation at cash flow dates only), the finite difference method 
can be faster by an order of magnitude than the Monte Carlo method to achieve 
the same accuracy in price. Finite difference method can be even more efficient 
in comparison with Monte Carlo in the case of local volatility model where Monte 
Carlo requires significantly larger number of time steps. In terms of robust and 
accurate estimation of Greeks, the advantage of the finite difference method will 
be even more pronounced. 



Keywords: Target Accumulation Redemption Note, option pricing, finite dif- 
ference, Monte Carlo 



1 Introduction 

Path dependent options have payoffs depending on the trajectory followed by one or 
more of the underlying processes. The most straightforward and easy to implement 
numerical solution for pricing path-dependent options is based on the Monte Carlo 
method. In the context of pricing path-dependent options by solving partial differ- 
ential equations (PDE), two additional challenges may merge due to the presence of 
path dependency. First, the dependency may introduce new dimensions to the partial 
differential equation. Second, it may cause the resulting equation much more difficult 
to solve because of the lack of diffusion in the additional dimensions. For some detailed 
discussions, see Tavella and Randall (2000), Zvan et al (1998) and Wilmott (2000b). 

The nature of the path-dependent option pricing problem largely depends on whether 
we have a continuous or discrete sampling for the path. In general, a continuous sam- 
pling model of path dependency introduces additional convection terms in PDE, while 
for a discretely sampled path-dependent option the convection terms are replaced by 
jump conditions. There are many successful attempts in pricing discretely sampled 
path dependent options by the PDE approach using lattice based method (e.g. bino- 
mial and trinomial trees used in Ritchken et al 1993, Hull and White 1993, Barraquand 
and Pudet 1996, Forsyth et al 2002), and similarly finite volume or finite element 
method (Forsyth et al 1999, Zvan et al 2001). Most of these studies consider Asian 
or lookback options. Typically, a linear interpolation is adapted in these methods in 
applying the jump conditions on the auxiliary variable (e.g. path average of the under- 
lying asset). The convergence study by Forsyth et al (2002) shows that it is possible 
for an algorithm based on lattice method to be non-convergent (or convergent to an 
incorrect answer) if the interpolation scheme is selected inappropriately. 

A Target Accumulation Redemption Note (TARN) provides a capped sum of pay- 
ments over a period with the possibility of early termination (knockout) determined 
by the target level imposed on the accumulated amount. A certain amount of payment 
(e.g. spot value minus the strike) is made on a series of cash flow dates (referred to 
as fixing dates) until the target level is breached. The payoff function of a TARN is 
path dependent in that the payment on a fixing date depends on the spot value of 
the asset as well as on the accumulated payment amount up to the fixing date. Typi- 
cally, commercial software solutions for pricing a TARN are based on the Monte Carlo 
method. This paper presents a finite difference scheme as an alternative to the Monte 
Carlo method to evaluate TARN. The focuses are on the step by step implementation 
of the finite difference scheme, which is not readily available in the literature, and on 
the comparison of performance of the proposed scheme relative to the Monte Carlo. 
We are not aware of any finite difference scheme published in the literature, although 



a general outline of PDE approach to pricing TARN can be found in Piterbarg (2004). 
Without losing generality, we assume the underlying asset is the foreign exchange 
(FX) rate. The definitions of TARN options with three different knockout types and 
some key notations are introduced in Section 2. Foreign exchange rate models are 
described in Section 3. Finite difference scheme for TARN is presented in Section 4 
and numerical results for both the finite difference and Monte Carlo methods are given 
in Section 5, before concluding by Section 6. 

2 TARN Payoff Definition 

There are different versions of TARN products used in FX trading. For simplicity, 
here we consider one specific form of TARN. The presented finite difference scheme 
can easily be adapted to other more general forms of TARN as discussed in Section 
14.2.31 Denote the FX rate at time t as S{t) and other notation as follows: to is today's 
date; K is the number of fixing dates (cash flow dates); ti,t2, . . . ,tK are fixing dates; 
X is strike; U is the target accrual level; S(ti), S(t2), ■ ■ ■ , Sitx) are FX rate values 
at fixing dates ti,t2, . . . ,tK', A{t) is accumulated amount at time t; and all amounts 
are per unit of notional foreign amount. On each fixing date tk, there is a cash flow 
payment 

Ck = PiS{tk) — X) X l/3x5'{tfc)>/3xX, (1) 

where /3 is a strategy on foreign currency {(3 = 1 corresponds to buy and /3 = — 1 
corresponds to sell), subject to the target level U is not breached by the accumulated 
amount A{tk). If the target level U is breached before or on the last fixing date, denote 
tj^ is the first fixing date when the target is breached, i.e. 

K = mm{k : A{tk) > U}, k = 1,2, . . . ,K . (2) 

Otherwise, set K = K. The actual payment on the fixing date tk < t^ can be written 
as 

Ck{S{tk), A(4„i)) = Ck^ (lA(t,_o+5fe<c/ + W^fe X lA(t,_o+c.>i/)' (3) 

and Cfc = for tk > tj^. Here, A{tk~i) is the accumulation amount immediately after 
the fixing date tfc_i, and Wk is the weight depending on the type of the knockout when 
the target level U breached. The accumulated amount A{t) is a piece- wise constant 
function A{t) = A{tk-i), t^-i ^ t < tk with 

A{tk) = A(t,„i) + Ck{S{tk), A{tk^,)). (4) 

There are three knockout types used in practice 



• Full gain - when the target is breached on a fixing date tk, the cash fiow payment 
on that date is allowed. This essentially permits the breach of the target once, 
and the total payment may exceed the target for full gain knockout. 

• No gain - when the target is breached, the entire payment on that date is disal- 
lowed. The total payment will never reach the target for no gain knockout. 

• Part gain - when the target is breached on a fixing date t^, part of the payment 
on that date is allowed, such that the target is met exactly. 

Formally, it can be represented by the following definition of the weight 

{1, if knockout type = full gain; 
0, if knockout type = no gain; (5) 

I3x(s(t )-x) ' ^^ knockout type = part gain. 
The present value (discounted value) of the TARN payoff in domestic currency for 
FX realization S = {S{ti), S(t2), ■ ■ ■ , Sitx)) is then 



fc=i 



Bd(tQ,tk) 



where [-Bd(to,^fc)]~^ is domestic discounting factor from the fixing date tk to to- 
Other forms of TARN used in trading include modifications of cash fiow payments 
([1]) and accumulated amount rule (jlj). In the present study, the cash fiow payment on 
each fixing date is the same as the increment in the accumulated amount, both are 
represented by Ck{S(tk),A{tk-i)). In other forms of TARN, the two quantities can 
differ, but this should cause no additional difficulties for the finite difference method 
presented here, as will be further discussed later in Section 14.2.31 

3 FX Model 

Under the standard no arbitrage option pricing methodology, today's fair price of 
TARN is calculated as the expectation of payoff (Q under the risk neutral process. 
Specifically, we consider the risk neutral process 

^ = (r, - rf)dt + adW,, (7) 

where r^ and rj are domestic and foreign local (instantaneous) interest rates, a is 
the local (instantaneous) volatility and Wt is the standard Brownian motion. The 
expectation can be calculated using Monte Carlo by simulating risk neutral process (ITj) 



many times and averaging the payoff realizations; or by solving corresponding PDE via 
the finite difference method. Here, the local interest rates can be constant or functions 
of time Td = rd{t), Vf = rf{t)] and volatility can be constant, function of time a = a{t) 
or function of time and FX rate a = a{S{t),t). The last case corresponds to the local 
volatility model that can be calibrated to match observed implied volatility surface; 
see e.g. Wilmott (2000a). 

4 Finite difference numerical scheme 

Let V{S,t,A) be the value of TARN for spot rate 5* and accumulated amount A at 
time t. Since the path-dependent quantity A is monitored discretely, there are no new 
diffusion terms and the standard option pricing PDE is still valid between fixing dates 

^ + -a\S, t)S'^ + (r,(t) - rf{t))S^ - r,(t)l^ = 0. (8) 

Typically, PDE solution for option pricing requires final conditions (the payoff) and 
boundary conditions (e.g. at zero or at a barrier). For discretely sampled path depen- 
dent options, additional jump conditions apply. Unlike in the case of path independent 
options where the payoff at expiry is known a priori and typically the final condition is 
of the Dirichlet type with the value of the payoff, in the case of the TARN option the 
final payoff is not known a priori. The expiry time is simply the last fixing time, and 
the final payoff depends on the path of the underlying up to the expiry time. Immedi- 
ately after the final payoff the option is worthless, and we can set the final condition 
to zero at T 

V{S,T,A) = 0, 

where T = t;^ is the last monitoring time. Applying a proper jump condition from 
T to T~ will give us a more informative final condition at T^, where T^ is the time 
infinitesimally before the last monitoring time tx = T. Unfortunately, any single 
solution of dS]) based on a given final condition at T" will not lead to the correct 
answer to the TARN option pricing, even if we know the final jump amount. We need 
multiple solutions to ([8]) with different final payoffs or jumps. Across any fixing date, 
there is a discontinuous but predictable jump in the accumulated amount. In such a 
case the no-arbitrage principle dictates that there must be a proper jump condition 
imposed on the path dependent option values. The jump value Ck given by (jlj) is the 
cash flow to the TARN owner, thus 

V{S, t,, A{t-)) = V{S, h, A{t,) + Ck{S, A{r,))) + Ck{S, A{t,)). (9) 

Finally, the PDE solution will give us the today's TARN price V{S{tQ),to,0). 



4.1 Jump condition application 

Let us introduce an auxiliary finite grid = Ai < y42 < ^3 ■ ■ ■ < Aj = U to track 
the accumulated amount A, where J is the total number of nodes in the accumulated 
amount coordinate. The upper limit U is needed because the accumulated amount 
cannot exceed the target U. For each Aj, we associate a continuous finite difference 
solution to the one- dimensional PDE (|8]). For finite difference solution, at every jump 
we let A to be one of the grid points Aj, 1 < j < J. Since A is always known at each 
jump to be one of the fixed nodal point values, there is no need to continuously track 
the actual evolution of the accumulated amount A during the entire finite difference 
solving process. 

Denote finite difference grid points in the S variable as 5*1, 5*2, ... , Sm, where M 
is the total number of nodes in the 5* coordinate. For any 5* = Sm, fn = 1, . . . , M, 
substituting v4(t^) with Aj, j = 1, . . . , J in (j9]) we get, 

V{Sm,t^,Aj) = V{Sm,h,A+) + Ck{Sm,A,), A] = A,+Ck{Sm,A,) (10) 

where t^ denotes the time infinitesimally before the monitoring time t^. In equation 
( ITOl) . we have let the accumulated amount before the k-th payment at t = t^ to be one 
of the grid point Aj. Equation (ITO|) describes a forward jump from t^ to tk- 

Because backward time marching is carried out for finite difference solution of PDE 
(jH]) associated with a fixed node point Aj, intuitively the jump should be applied 
backwards from tk to t^ . That is, in finite difference solution the value of A at tk is 
known to be one of the grid point Aj, and after a backward jump from tk to t^ the 
value of A changes from Aj to Aj . This backward jump can be expressed as 

V{Sm, tl, A-j) = V{Sm, tk, Aj) + Ck{Sm, A-j), A-j = Aj - C,(5^, Aj) (11) 

In both ( ITOl) and (iTTl) . Ck{S,A) is calculated according to (jS]). Figure 1 illustrates 
the application of jump condition (TTTl) . 



4.2 Tracking finite diff^erence solutions 

The idea is tracking J finite difference solutions corresponding to the J grid points for 
the auxiliary variable, the accumulated amount. For each fixed accumulated amount 
Aj, ^ <i j <i J 1 we start solving PDE by the finite difference scheme with the final 
condition V{S,T,A) = and a final jump condition from tx = T to tj. = T~. The 
implementation of this idea is not straightforward, because at each sampling time, 
the jump condition ( TTTl) has to be applied and the accumulated amount after each 
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t = t 



t = t 




Figure 1: Illustration of jump conditions applied to finite difference grids. 



jump changes accordingly and falls off the grid points Aj, I < j < J. Not only the 
accumulated amount changes, the amount of changes differ for different grid points in 
the underlying space. As shown in ( iTTl) . because Aj is not a constant, the solution 
V{Sm,t'i^,A~) obtained after the jump does not correspond to any grid point in the 
auxihary variable space. Worse still, the set of values V{Sm,t~,Aj), m = 1, . . . ,M 
does not correspond to any continuous finite difference solution of the one-dimensional 
pde - it does not satisfy the PDE because the value Aj is scattered all over the place, 
not associated with any unique value. This is because for the option value to satisfy 
the one-dimensional pde, it requires a unique accumulated amount at any time - for 
consistency one cannot have different accumulated amounts at the same time. 

For the M grid points in S space, V{Sm, t'^, ^J), ""^ = I, ■ ■ ■ ,M correspond to M 
different scenarios of payoffs. Before the jump, the M values V{Sm, tk, Aj) are related to 
each other through the PDE, because they are all associated to the same accumulated 
amount Aj. The connection between the M values is broken after the jump. 

4.2.1 Reversal of the jump direction 

Intuitively, jump conditions should be applied through fITT]) . That is, as the backward 
marching is performed for each of the J solutions corresponding to Aj, 1 < j < J, 
at any crossing of sampling time from tk to t^, Aj jumps to A~ and the solution 
V{Sm, tk, Aj) jumps to V{Sm, t^, A~) according to (TTTI) . We can then interpolate from 
V{Sm,tk^ ^]) ^'^ obtain V(S'm,t^,^j) and continue time marching backwards until 
next sampling date. 

Unfortunately, the intuitive application of jump conditions as described in the above 
paragraph is problematic in two important ways. First, from flTTj) . it is possible to get 
a negative value for A^ , which is invalid (meaningless) and out of the range of the 
auxiliary variable space; Second, because Q < Aj <U and flTTj) for Aj is a decreasing 
function, the target U will never be exceeded by any of the jump according to ( ITTj) . thus 
there is no way to apply the different knockout conditions for the different knockout 
types as specified in ([5]). In other words, applying jump condition using (ITT]) cannot get 
the correct answers to any of the knockout types of TARN. Essentially, applying ( ITTj) 
artificially restricts the boundary for the auxiliary variable to be within the target, 
instead of letting the underlying process breach the target. Another minor issue is 
that the jump condition f lTTj) is implicit in A^, i.e. strictly speaking the jump amount 
Cfc(5'm, A~) is not known before the backward jump. 

The remedy to the above problems is actually quite simple - we reverse the direction 
of the jump. Jump condition ([9]) is true for any values of the auxiliary variable in the 
range < A < [/, i.e. we do not have to use Aj, a grid point, on the right hand side 



8 



as in (TTTI) . Instead, we could use (fTOl) to have the value of A at t^ to be one of the 
gird point Aj. Now (TTOj) is explicit in Aj. Since grid point Aj satisfies < Aj < U 
and A^ > Aj, after the jump from Aj to A^ , A^ will never be negative. What is 
more, A^ value may now exceeds the target U, allowing the knockout conditions to be 
imposed. The knockout condition is implied in the calculation of Ck{Sm,Aj) in fITOl) . 
using equation ([3]). Specifically, we have 

Ck{S^, Aj) = CkX (Ia^^c.ku + WkX ^A,+c,>u)^ (12) 

Equation ( ITOll gives the desired solutions V{Sm,t'^ ,Aj) at the gird points Aj, j = 
1, . . . , J, given V{Sm, tk, A^). Because we only have solution V{Sm, tk, Aj) upon march- 
ing to time tk, we need performing interpolation from V{Sm, tk, Aj) to obtain V{Sm, tk, A^ 
for all m = 1,2, . . . ,M and j = 1,2, . . . , J. 

4.2.2 Cubic Spline interpolation 

For a fixed grid point in spot S = Sm, there are J values after the jumps corresponding 
to the J solutions associated with each of the J gird points in the auxiliary variable. 
These values are given by V{Sm,tk,Aj), j = 1, . . . ,J. We need to extract J values 
V{Sm,tk, A^) from V{Sm,tk,Aj), j = 1, . . . , J by interpolating with respect to Aj. For 
a given 5*^ this is a one-dimensional interpolation in the accumulated amount space. 

As shown in a convergence study by Forsyth et al (2002), it is possible for a numerical 
algorithm for discretely sampled path-dependent option pricing to be non-convergent 
(or convergent to an incorrect answer) if the interpolation scheme is selected inappro- 
priately. All the previous studies of numerical PDE solution for path dependent (Asian 
or lookback options) used either a linear or a quadratic interpolation in applying the 
jump conditions. In our experience a better choice is the cubic spline interpolation 
(Press et al 1992). This procedure assumes the J values, V{Sm,tk,Aj),j = 1, . . . , J, 
form a smooth function in the auxiliary variable space and the cubic spline interpola- 
tion has a much higher order of accuracy than linear or quadratic interpolation. The 
error of cubic spline is 0{h^) where h is the size for the spacing of the interpolating 
variable, assuming a uniform spacing. In our case h = 6 A = U/{J — 1). Natural 
boundary conditions are imposed at the two ends Aq = and Aj = U, i.e. we assume 
zero second derivative of the spline function at the two ends. For each fixed spot Sm, 
a single tri-diagonal system of equations is solved once for obtaining all the J values 
V{Sm,tk,A+), j = l,...,J. 

If we perform the above interpolation for all the M grid points in spot 5* and apply 
jump condition (ITOl) . we will have M x J new values 1^(5*^,^^,^^), rn = 1,...,M, 
j = 1, . . . , J. For a fixed j, the M new values V{Sm, t'^ , Aj) correspond to the PDE 
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solution associated with grid point Aj. Given V{Sm,t'j^,Aj) for each fixed Aj, we 
can now continue time marching backwards until the next sampling time. The whole 
algorithm can be summarized as follows. 

1. Apply zero final condition a.t T = tx for all the J solutions to equation (jH]) 
corresponding to Aj, I < j < J. 



2. Apply the jump condition ( jTOl) to obtain A^ for each of the J solutions at each 
of the M grid points in spot, beginning with k = K [t^ = T) for the first jump. 

3. Perform cubic spline interpolation from points V{Smi tk, Aj),j = 1, . . . , J to new 
points V{Sm, tk, A'j') by forming a smooth function from the J values V{Sm, tk, Aj) 
with each spot grid point 5*^. 

4. Apply the jump condition (ITU1) . i.e. calculate V{Sm,t^ ,Aj) from V{Sm,tk,A^). 
The knockout condition (a boundary condition in variable A) is implied by the 
calculation of Ck{Sm,Aj) using flT2|) . 

5. Perform the finite difference time marching backwards for each of the J solutions 
V{S,t^ ,Aj), j = 1, . . . , J, corresponding to the J grid points in the auxihary 
variable, until a sampling time is encountered. This gives solution V{S, tk-i, Aj). 

6. Repeat steps 2 to 5 until k = 1. 

7. Take the single solution V{S,ti,0) to do final time marching until t = to, and 
take V{S, to, 0) as the final solution of the TARN option. 

As indicated in Step 1, at the final fixing time T = t^, the zero- value final condition 
is applied at T and the following jump condition is applied before taking any pde 
solving steps 

V{Sm, T-, A,) = + Ck{S^, A,). (13) 

For each set of V{Sm,T~ ,Aj) with fixed j, we begin tracking a finite difference 
solution through backward time marching. 

In step 7, only a single solution is needed between the first sampling time t^and the 
spot date - there is no more need to track all J solutions, since there are no more jump 
conditions to be applied. For good accuracy, we require that the current spot value 
^(to) be one of the grid point in S = {Si, 5*2, ... , Sm)- 
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4.2.3 Extension to other TARN products 

The payoff structure with the three k;nock;out types of TARNs considered in this study 
is typical in FX trading, but there are other types with different payoff structure and 
knockout type. The extension of the present FD method to other TARN types is 
straightforward. For example, suppose there are extra payments CI at each fixing date 
tk and this extra payment does not count in the knockout condition ([3]) but wiii also 
get knockout by the same knockout condition ([3]), i.e. 

Cl{S,A{t,.,)) = CI X (1^. _,Hc;<c/ + Wk X Imu_.,-,K>u)^ (14) 



K 



{Cu (5(tfc), A{tk-i)) + CI (5(4), A{tk-i))) 



p^Q^^y y^kyoy^kj,^y'^k-w-r^,yoy.,j,^y.,-,jjj ^ A(to) = 0, (15) 

where C^ is the extra payment when the target is not breached. In this case, the only 
change in the finite difference scheme is to replace the price jump condition (TTOl) with 
a new condition 

V^(5„,t„ A,) = F(5^,4, A+) + Cu{S^,A^) + C*(5„,A,). (16) 

There is no any other change required in dealing with the auxiliary variable Aj, 
since the extra payment does not contribute to the monitored accumulated amount A 
and the knockout condition remains the same. 

4.3 Boundary condition 

Typically a finite difference solution is sought within a rectangular domain (0 < t < 
T, < Srain ^ S < ^max), whcre both Srain and ^max are chosen to be sufficiently 
far away from the spot price of the underlying asset, e.g. three standard deviations 
from the spot. To insure an unique solution, boundary conditions are required at S'mm 
and S'max- There are different ways of imposing proper boundary conditions that are 
numerically equivalent. A rather general and robust boundary condition at both 5*111111 

and Omax IS 

(,*-'min) t-j — U, ^ (,»-'niax) t-j — U, 



which is particularly useful because it is independent of the contract being valued, 
provided the option has a payoff that is at most linear in the underlying for small and 
large values of S (almost all common contracts have this property). Other boundary 
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conditions work equally well. For example, for a call option, the following boundary 
condition can be applied 

dV 

V (Omiri) t) = U, -;r— (Omax? t) = 1; 



dS 



and for a put option we have 



dV 

^ ('-'max; t) = U, -j-^ ('-'mm) t) = ~1- 

Some detailed discussions on various suitable boundary conditions can be found in 
Wilmott (2000b). 

4.4 Log-transform 

It is a common practice to re-write equation ([H]) in terms of x = In(S') before finite 
difference discretization: 

dV 1 ^dV^ dV 

where v = rait) — rjif) — (t^/2. Equation (1T7|1 is slightly simpler than ([S]), i.e. if 
volatility and interest rates are constant, then coefficients of all derivatives in (IT7|) are 
all constant. 



4.5 Discretization for uniform grid 

Unlike barrier options, pricing the discretely monitored TARN option can always rely 
on uniform grids. This is because there are at most two critical points to be 'pinned' 
to grid points - the spot and the strike, provided we make the far boundaries flexible. 
Since the only requirement for far boundaries is that they are sufficiently far from spot, 
these boundaries can certainly be extended a bit further to accommodate uniform grids 
with the two critical points (spot and strike) pre-determined. When the spot and the 
strike are almost the same, uniform grids tied to both the spot and the strike may have 
too large a number of nodes, in this case we chose to tie the strike only, and perform 
a one-off final interpolation to obtain the price corresponding to the spot. 

Denote the option price at time step n and grid point Si as V^, n = 0,1,2, ... ,N. 
For a uniform grid, 6xi = Xi — Xi^i = 6x is a. constant, and we obtain the following 
finite difference approximation with second order accuracy 

^{x.,tn)= '\^'~' +0{6x'), (18) 
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a2T/ yn _ csyn , yn 



The 6'— scheme 

Define the following differential operator -F(V, x, a, z/, r^) 



1 (9\/^ (9\/ 
F(V^, X, a, ^, r,) = -a^^^ + z.— - r,y (20) 



and the associated finite difference operator F" 



-o-^(xi,t„)^p^(xi,t„) + z/(xi,t„)--(xi,t„) -rd(t„)\/.", (21) 



1 2. s^. ^^ 

where the first and second derivatives are approximated by finite difference as discussed 
above. Then the 6*— scheme can be expressed as 



yn+l _yn 

^-— '- + eFp^^ + (1 - 9)Fp = 0, (22) 

where < ^ < 1. Special values of 6* = 0, 6* = 0.5 and 9 = 1 correspond to fully 
explicit, Crank- Nicholson and fully implicit scheme, respectively. 



5 Numerical examples 

Comparison of the finite difference and Monte Carlo methods is performed in the case 
of basic model with constant volatility. In this case, the number of time steps for Monte 
Carlo simulated paths is the same as the number of fixing dates. In the case of basic 
or term structure models, simulations between fixing dates are not required because 
transition density between fixing dates is known in closed form (it is just a lognormal 
density). For local volatility model, simulations between fixing dates are required that 
will increase computations proportionally to the number of time steps. 

In the examples we consider all three types of knockout as described in Section 
2, each knockout type has four cases with four different targets, so the total number 
of numerical examples is 12. The other inputs common to all the examples are spot 
^(O) = 1.05, strike X = 1.0, volatility a = 0.2, interest rates r^ = rf = 0, fixing dates 
are every 30 days and we assume 20 fixing dates. 

Results are summarised in Table [H As shown in Table [H the computing time for 
Monte Carlo estimates based on Nsim = 200, 000 simulated paths is very close to that 
for the finite difference method with mesh 500 x 100 x 500 (500 points for spot, 100 
points for accumulated amount and 500 steps for time). 

In Table [Tj the Monte Carlo standard error is compared with estimated relative 
error of the finite difference solution. Ideally, relative error should be computed as the 
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relative difference between numerical solution and the exact solution, for both Monte 
Carlo and finite difference. Unfortunately in the case of TARN options, closed form 
solution cannot be found except limiting cases of one fixing date or very large target 
level. Nevertheless, the standard error in Monte Carlo and the estimated relative error 
in finite difference are both very good approximate to the exact relative error. In the 
case of finite difference, we estimate the relative error by using solution of the refined 
grids in spot, accumulated amount spaces as well as in time. Specifically, we double 
the number of grid cells in all three dimensions for the refined calculation, i.e. using 
grids 1000 x 200 x 1000 for spot, accumulated amount and time, and use this refined 
solution in place of the exact solution in estimating the relative error. As shown in 
the Appendix, because the 6*— scheme is second order in accuracy in both spot space 
and time, and the cubic spline interpolation in the accumulated amount is of the order 
0(/i^), using the solution of the refined grids in estimating the true relative error of 
the coarser grids is valid and well justified. 



target 


MC 


FD 


difF% 


stderr MC % 


MC sec 


err FD % 


FD sec 


No gain 


0.3 


0.1955 


0.1955 


0.0000% 


0.10% 


1.31 


0.045% 


1.12 


0.5 


0.3288 


0.3286 


0.0609% 


0.10% 


1.32 


0.001% 


1.13 


0.7 


0.4507 


0.4505 


0.0443% 


0.10% 


1.32 


-0.018% 


1.13 


0.9 


0.5633 


0.5633 


0.0000% 


0.10% 


1.32 


0.015% 


1.14 


Part gain 


0.3 


0.2446 


0.2445 


0.041% 


0.08% 


1.32 


0.016% 


1.12 


0.5 


0.3819 


0.3818 


0.0262% 


0.09% 


1.33 


0.005% 


1.13 


0.7 


0.5063 


0.5061 


0.0395% 


0.10% 


1.32 


0.038% 


1.13 


0.9 


0.6203 


0.6200 


0.0483% 


0.10% 


1.32 


0.010% 


1.13 


Full gain 


0.3 


0.2979 


0.2978 


0.0336% 


0.08% 


1.32 


0.039% 


1.12 


0.5 


0.4389 


0.4386 


0.0684% 


0.09% 


1.33 


0.001% 


1.12 


0.7 


0.5646 


0.5644 


0.0354% 


0.10% 


1.33 


0.015% 


1.13 


0.9 


0.6792 


0.6790 


0.0295% 


0.10% 


1.32 


0.012% 


1.13 



Table 1: Finite Difference (FD) vs Monte Carlo (MC) results for TARN price. The 
notional amount is one unit of foreign currency. The column "diff %" shows the relative 
difference between results of FD and MC. The computing time is for desktop with Intel 
Core 15-2400 @3.10GHz and 4 Gb RAM. 

As show in Table [H, the accuracy of finite difference solution is significantly better 
than that of the Monte Carlo in all the 12 test cases. On average, the Monte Carlo 
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standard error is about 0.1%, while the finite difference relative error is about 0.02%. 
That is, Monte Carlo relative error is five times as large as finite difference relative 
error, thus on average Monte Carlo computing time should increase by the factor of 25 
to achieve the same accuracy as finite difference because Monte Carlo standard error 



is proportional to l/^JNsim- Note that quoted Monte Carlo relative error is computed 
from the standard error of the estimate, i.e. it should be at least doubled for a more 
realistic error estimate. To improve the accuracy of Monte Carlo estimates, in our 
numerical example, we use the sum of vanilla options with maturities at the fixing 
dates as a control variate error reduction technique. Monte Carlo efficiency can also be 
improved by the use of other error reduction techniques such as importance sampling 
described in Piterbarg (2004) but it might be difficult to implement this for more 
general models such as local volatility model and we did not pursue this further. 

These numerical results clearly demonstrate that the use of finite difference will 
be even more beneficial (in terms of accuracy) in the case of local volatility model 
where Monte Carlo method will require simulations for extra time slices between fixing 
dates. We expect that the impact in efficiency will be more pronounced in calculation 
of Greeks where even small error in price such as 0.1% may lead to 10-100% error in 
second derivatives (e.g. Gamma or Vanna). 

6 Conclusions 

We have implemented a finite difference scheme for evaluating TARN options. Numer- 
ical results show that finite difference scheme is more efficient in pricing TARN than 
the Monte Carlo counterpart, even for basic models where the volatility is constant or 
piecewise constant between fixing dates. For a surface model, the computing time in 
the Monte Carlo method will increase in proportion to the number of time steps in 
the surface model, while the finite difference scheme presented here remain essentially 
the same in terms of computing time. In the numerical examples only price was con- 
sidered. It is expected that if the Greeks are considered in the comparison between 
finite difference and Monte Carlo, the advantage of finite difference will be much more 
significant. Even a small error in price such as 0.1% may lead to a large error 10-100% 
in second derivatives (e.g. Gamma or Vanna). Thus pricing TARN and its Greeks by 
the proposed finite difference scheme provides significant practical advantage over the 
commonly used Monte Carlo method. 

We have given very detailed descriptions of the numerical steps required in the finite 
difference scheme, so that readers can easily follow the procedures to implement their 
own, and re-produce the result if desired. The TARN structure considered in this 
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study is simple. However, implementation of the finite difference method can be easily 
extended to a more generalised accumulation rule and TARN parameters varying across 
fixing dates. 

7 Appendix: Estimation of Numerical Error 

Denote as follows: V is the exact solution, V is the numerical solution of the coarser 
grids (e.g. 500 x 100 x 500 for spot, accumulated amount and time), V* is the numerical 
solution of the refined grids doubled in each direction (e.g. 1000 x 200 x 1000) and 
S = V — V and 6* = V* — Vare the absolute numerical errors of the two grids, 
respectively. Then the relative difference between the numerical solutions of the coarse 
girds and refined grids is 



V -V* 
V* 


= 


6-6* 

V* 



and the true relative difference between the numerical solution of the coarse girds and 
the exact solution is 



V -V 



V 



It is easy to show e is a very good approximation of e. Due to the second order 
accuracy in both time and space, and fourth order accuracy in the accumulated amount 
cubic spline interpolation, 6* can be estimated as 6* = ±2~^6. Thus 



6-6* 



V* 



5^2-^5 



V + 2-^6 



where the last approximation sign is due to 



6_ 

V 

256=Fl 



256^1 



256 + 6/V 



1, assuming \6\ « V. Depend- 



256+<5/y 

ing on the signs of the absolute errors, e could be slightly overestimating or slightly 
underestimating the true relative error e. Thus using relative error between solutions 
of coarser grids and the refined grids (with number of grids in all dimensions doubled) 
as an estimate of the true relative error is well justified for a numerical scheme with 
second order accuracy. 



16 



